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Comment on "Twisted Protein Aggre- 
gates and Disease: The Stability of Sickle 
Hemoglobin Fibers" 

In a recent paper pQ Turner et al constructed the free 
energy per unit volume, G, needed to create a fiber bun- 
dle, where G = F— using continuum elasticity theory. 
Here F is the distortion free energy per unit volume of 
a bundle of radius R and pitch length A and ip is a pos- 
itive Lagrange multiplier that controls the pitch length. 
From G they predicted the physical properties of the fiber 
bundle, such as the equilibrium (metastable) bundle ra- 
dius R c , where in the latter case they minimized G with 
respect to R. However, we believe their analysis is incor- 
rect for two reasons, the first being the use of the free 
energy density, G, rather than the total free energy O, to 
determine R c . The second is their omission of the bind- 
ing energy between fibers, which in classical nucleation 
theory of spherical droplets corresponds to the driving 
force for nucleation. We present a corrected version of 
their analysis below. Our approach is the same as that 
of Grason and Bruinsma [2], who determined the critical 
bundle size for aggregates of filamentous actin. 

According to classical homogeneous nucleation theory 
[3 0] , the critical " droplet" size corresponds to the mini- 
mum of the total free energy f2 = R 2 LG, which is signif- 
icantly different from the energy density G. A simple ex- 
ample is the nucleation of a spherical droplet [21 0j ■ The 
analogous argument for the heterogeneous nucleation of 
the fiber bundle involves calculating the total free energy 
required to create this bundle from an aggregate of fibers 
of (fixed) length L. 

The energy of a twisted fiber as a function of pitch A 
and radius R, includes the contributions from the surface 
tension, extension or compression, bending, twisting and 
binding: 
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where E is the extensional modulus, a the radius of a 
protofilament and L the fiber length, ip is related to the 
twisting stiffness [IJ[2]. Equation[I]contains an additional 
term — R 2 e due to the aggregation energy 00] between 
fibers that is not present in Turner et al. pQ. The equi- 
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Using experimental values for HbS of a = Anm, 
E = 51MPa, V = 3.5 x l^JmT 2 p][5], we find that 
n(R;A*) has just a single peak for e > 190Jm~ 3 (Fig. 
[l}. R = and R — >• oo correspond to the phases of the 
dispersed protofilaments and stable crystal structures, re- 
spectively. As e decreases below this, a local minimum 
develops in f2(i?;A*) whose position depends on e and 
7. The minimum critical bundle size R r occurs under 


s = 107m 3 ,y = 5.4/A/m -2 
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FIG. 1: (Color online) Plot of the free energy f2(7?;A*) 
per unit length as a function of the fiber radius, R, us- 
ing the experimental values for HbS given in the text for 
e > 190 Jm" 3 (black), e = 38 Jm" 3 (red), and e = 10 Jm" 3 
(green). One local minimum occurs at R — llnm which cor- 
responds to the (metastable) equilibrium radius of HbS when 
e = 38 Jm" 3 , 7 = 5.8/iJm" 2 . 


Combining the estimate e s» 38 Jm~ 3 for HbS [S], this 
yields a value of R c = llnm and 7 = 5.8/iJm~ 2 ( Fig. 
nl), which are consistent with experimental observations 
[Tj [S] . A further reduction in e leads to a decreasing value 
of fl(R c ;A*) (Fig. [I]). We also note that the torsional 
rigidity obtained by Turner et al is the same in our cal- 
culation and in agreement with experimental values. Fi- 
nally, there always is an energy barrier for the transition 
from dispersed protofilaments to the metastable bundle 
phase, which is incorrectly predicted as a spontaneous 
process in reference [T]. 
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